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We study the renormalization group flow of the Lagrangian for statistical and quantum systems 
by representing their path integral in terms of a tensor network. Using a tensor-entanglement- 
filtering renormalization (TEFR) approach that removes local entanglement and coarse grain the 
lattice, we show that the resulting renormalization flow of the tensors in the tensor network has a 
nice fixed-point structure. The isolated fixed-point tensors Ti nv plus the symmetry group G S ym of 
the tensors (ie the symmetry group of the Lagrangian) characterize various phases of the system. 
Such a characterization can describe both the symmetry breaking phases and topological phases, as 
illustrated by 2D statistical Ising model, 2D statistical loop gas model, and 1+1D quantum spin- 
1/2 and spin-1 models. In particular, using such a (G S ym, Tinv) characterization, we show that the 
Haldane phase for a spin-1 chain is a phase protected by the time-reversal, parity, and translation 
symmetries. Thus the Haldane phase is a symmetry protected topological phase. The (G sym ,Ti nv ) 
characterization is more general than the characterizations based on the boundary spins and string 
order parameters. The tensor renormalization approach also allows us to study continuous phase 
transitions between symmetry breaking phases and/or topological phases. The scaling dimensions 
and the central charges for the critical points that describe those continuous phase transitions can 
be calculated from the fixed-point tensors at those critical points. 
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I. INTRODUCTION 

For a long time, people believe that all phases and all 
continuous phase transitions are described by symmetry 
breaking. 1 Using the associated concepts, such as order 
parameter, long range correlation and group theory, a 
comprehensive theory for symmetry breaking phases and 
their transitions is developed. However, nature's richness 
never stops to surprise us. In a study of chiral spin liq- 
uid in high temperature superconductor, 2 ' 3 we realized 
that there exist a large class of phases beyond symme- 
try breaking. 4 Since the new phases can be characterized 
by topology-dependent topologically stable ground state 
degeneracies, the new type of orders in those phases are 
named topological orders. 5 . Later, it was realized that 
fractional quantum Hall states 6 are actually topologically 
ordered states. 7 

Recently it was shown that topological order is actu- 
ally the pattern of long range entanglement in the quan- 
tum ground state. 8-10 The possible patterns of long range 
entanglement are extremely rich, 11-13 which indicate the 
richness of new states of matter beyond the symmetry 
breaking paradigm. Such rich new states of matter de- 
serve a comprehensive theory. 14 ' 15 

But how to develop a comprehensive theory of topo- 
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logical order? First we need a quantitative and con- 
crete description of topological order. The topological 
ground state degeneracy 4,7 is not enough since it does not 
completely characterize topological order. However, the 
non-Abelian Berry's phases from the degenerate ground 
states and the resulting representation of modular group 
provide a more complete characterization of topological 
order. 5 It may even completely characterize topological 
order. There are other ways to (partially) characterize 
topological order, such as through the fractional statistics 
of the quasiparticles 16 ' 17 or through the string operators 
that satisfy the "zero-law" . 18 Those "zero- law" string op- 
erators appear to be related to another characterization 
of topological order through a gauge-like-sysmetry de- 
fined on a lower dimensional subspace of the space. 19 ' 20 

But so far, we have not been able to use those charac- 
terizations to develop a theory of topological order. So 
the key to develop a comprehensive theory of topological 
order is to find a quantitative and easy-to-use description 
of topological order. For fractional quantum Hall states, 
the pattern of zeros appears to provide a quantitative 
and easy-to-use description of topological order. 21-23 For 
more general systems, we may use the renormalization 
group (RG) approach to find a quantitative and easy- 
to-use description. We will see that the new description 
based on the RG approach naturally leads to a theory of 
topological order and phase transitions between topolog- 
ically ordered states. 

Renormalization group (RG) approach 24,25 is a pow- 
erful method to study phases and phase transitions for 
both quantum and statistical systems. For example, for a 
quantum system, one starts with its Lagrangian and path 
integral that describe the system at a certain cut-off scale. 
Then one integrates out the short distance fluctuations 
in the path integral to obtain a renormalized effective 
Lagrangian at a longer effective cut-off scale. Repeating 
such a procedure, the renormalized effective Lagrangian 
will flow to an isolated fixed-point Lagrangian that rep- 
resents a phase of the system. As we adjust coupling 
constants in the original Lagrangian, we may cause the 
renormalization flow to switch from one fixed point to 
another fixed point, which will represent a phase tran- 
sition. Since the fixed-point Lagrangian is unique for 
each phase and the phase transition is associated with a 
change in fixed-point Lagrangian, it is natural to use the 
fixed-point Lagrangian to describe various phases. 

But can we use the fixed-point Lagrangian to quanti- 
tatively describe topological order? At first sight, the an- 
swer appears to be no. When a quantum spin system (or 
a qbit system) has a topologically ordered ground state, 
its low energy effective Lagrangian {it the fixed-point La- 
grangian) can be a Z2 gauge theory, 26-30 a f/(l) gauge 
theory, 29,31-33 a SU(2) gauge theory, 29 a QED theory 
(with U{1) gauge field and massless fermion field), 29,34,35 
a QCD theory (with non-Abelian gauge fields and mass- 
less fermion fields), 29,34,36,37 a Chern-Simons theory, 3,38 
etc . In the standard implication of RG approach (within 
the frame work of field theory), a bosonic Lagrangian 



that describes a quantum spin system cannot flow to all 
those very different possible fixed-point Lagrangian with 
very different fields (such as gauge fields and anticom- 
muting fermionic fields). 

Since lattice spin systems do have those very different 
topologically ordered phases, we see that the standard 
implementation of the RG approach is incomplete and 
needed to be generalized to cover those existing (and ex- 
citing) renormalization flows that end at very different 
fixed-point Lagrangian. In this paper, we will describe 
one such generalization - tensor network renormalization 
approach - that should allow us to study any phases 
(both symmetry breaking and topological phases) of any 
systems (with a finite cut-off). The resulting fixed-point 
Lagrangian from the tensor network renormalization ap- 
proach give us a quantitative description of topological 
orders (and the symmetry breaking orders). 

The tensor network renormalization approach is based 
on an observation that the space-time path integral of a 
quantum spin system or the partition function of a sta- 
tistical system on lattice can be represented by a tensor 
trace over a tensor network after we discretize the time 
and the fields: 

f Dip(x,t)c-f c = J2U T ^^ =™®i T (!) 

{00 

where i label the discretized space-time points and <f>i is 
the discretized field (for detail, see (2)). The different 
choices of tensor T^ >i ^.^ k ^ >l in the tensor network corre- 
spond to different Lagrangians/partition- functions. We 
like to point out that in addition to use it to describe path 
integral or partition function, tensor network can also 
be used to describe many-body wave functions. 39 Based 
on such a connection, many variational approaches were 
developed. 39-45 

Since every physical system is described by path inte- 
gral, so if we can calculate the tensor trace, we can cal- 
culate anything. Unfortunately, calculating the tensor- 
trace tTr is an NP hard problem in 1+1D and higher 
dimensions. 46 In Ref. 47, Levin and Nave introduced a 
tensor renormalization group procedure where a coarse 
graining of the tensor network is performed to produce a 
new tensor on the coarse grained network. This leads to 
an efficient and accurate calculation path integral (or par- 
tition function) and other measurable quantities. Such a 
method opens a new direction in numerical calculation 
of many-body systems. 

The above tensor renormalization group procedure also 
gives rise to a renormalization flow of the tensor. How- 
ever, the fixed-point tensor obtained from such a flow 
is not isolated: even starting with different Lagrangian 
{it the tensors) that describe the same phase, the renor- 
malization flow of the those Lagrangian will end up at 
different fixed-point tensors. This means that the renor- 
malization procedure does not filter out (or integrate out) 
all the short distance fluctuations so that the fixed-point 
tensors still contain short distance non-universal informa- 
tion. In this paper, we will introduce an improved ten- 
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sor renormalization group procedure by adding a filter- 
ing operation that remove the local entanglements. We 
show that this additional entanglement filtering proce- 
dure fixes the above problem, at least for the cases stud- 
ied in this paper. 

To test the new tensor network renormalization 
method, we study the fixed-point tensors that correspond 
to the two phases of 2D statistical Ising model which has 
a Z2 spin flip symmetry. We find that the fixed-point ten- 
sor for the disordered phase is a trivial tensor T TRI where 
the dimensions of all the tensor indices are equal to one 
(we will refer such a trivial tensor as a dimension-one ten- 
sor which is represented by a number). While the fixed- 
point tensor for the symmetry breaking (spin ordered) 
phase is a dimension-two tensor T Z2 which is a direct 
sum of two dimension-one tensors: T Z2 = T™ T™. 
We believe that this is a general feature for symmetry 
breaking phase: a symmetry breaking phase with n de- 
generate states is represented by fixed-point tensors that 
is a direct sum of n dimension-one trivial tensors. 

Since the Z 2 spin flip symmetry is crucial for the ex- 
istence and the distinction of the disordered phase and 
the spin ordered phase, we should use a pair, the sym- 
metry group G sym = Z 2 and the fixed-point tensors 
Tinv = T™, T z ' 2 , to characterize the two phases. So 
the disordered phase is characterized by (Z 2 ,T TRI ) and 
the spin ordered phase is characterized by (Z 2 , T Z2 ) . It 
turns out that the characterization based on the symme- 
try group and the fixed-point tensor, (G sym , T; nv ), is very 
general. It can describe both symmetry breaking phases 
and topological phases. This is one of the main results 
of this paper. 

Even at the critical point between the symmetry break- 
ing phases and/or topological phases, the tensor renor- 
malization approach produces fixed-point tensors (which 
now have an infinite dimension). From those fixed-point 
tensors, or more precisely, from a finite-dimension ap- 
proximation of those fixed-point tensors, we can calculate 
the critical properties of the critical points with an ac- 
curacy of 0.1%. Using such an method, we confirm that 
the phase transition between the disordered and ordered 
phases of the Ising model is a central charge c = 1/2 
critical point. 

We also study a 2D statistical loop gas model. We 
find that for the small-loop phase, the fixed-point ten- 
sor is the dimension-one tensor T™, while for the large 
loop phase, the fixed-point tensor T LL is a tensor with 
dimension 2 (see (24)). However, unlike the symmetry 
breaking phase, the fixed-point tensor for the large loop 
phase is not a direct sum of two dimension-one tensors: 

T LL T TRI e T TKl Thig ig consistent w j th the f act 

that, as a statistical system, the small- loop phase and the 
large loop phase has the same symmetry and the phase 
transition is not a symmetry breaking transition. The ex- 
ample demonstrates that fixed-point tensors can be used 
to study phases and phase transitions beyond symmetry 
breaking. The phase transition between the two phases 
is again a central charge c = 1/2 critical point. 



Next, we study a spin-1/2 model H = J2i(~ a f a i+i ~ 
•/ofof+i), with Z% x Zf symmetry generated by ia x 
and ier z . The spin-1/2 model has two different Z 2 sym- 
metry breaking phases described by, in terms of the 
new notation, (Zf x Zf,T Z2 ) and (Zf x Zf,T LL ). The 
(Zf x Zf , T Z2 ) phase breaks the Zf symmetry while the 
(Zf x Zf , T LL ) phase breaks the Zf symmetry. The phase 
transition between the two Z 2 symmetry breaking phases 
is a continuous phase transition that is beyond the Lan- 
dau symmetry breaking paradigm. We find that the ten- 
sor network renormalization approach can describe such 
a continuous phase transition. The phase transition be- 
tween the two phases is a central charge c = 1 critical 
point. 

Last we study the T = phases of quantum spin-1 
model in 1+1D H = -^,[5, • S i+1 + U(Sf) 2 + BSf], 
and found that the model has four different phases for 
different (U, B): a trivial phase, two Z 2 symmetry break- 
ing phases, and the Haldanc phase. 48-50 We find that the 
Haldane phase is described by a dimension-four fixed- 
point tensor T Haldanc (see (31)). T Haldanc is not a direct 
sum of dimension-one tensors. From the fixed-point ten- 
sor yHaidano^ wc confirm that the Haldane phase has a 
finite energy gap and the ground state is not degenerate. 
Thus, despite the fixed-point tensor to be a dimension- 
four tensor, the Haldane phase, like the trivial phase, 
does not break any symmetry (as expected). 

The TEFR approach allows us to study the stabil- 
ity of the Haldane phase (ie the distinction between the 
Haldane phase and the trivial phase) under a very gen- 
eral setting. Our calculation shows that the fixed-point 
tensor J lHaldane ; and hence, the Haldane phase is robust 
against any perturbations that have time-reversal, parity, 
and translation symmetries. 73 So the Haldane phase is a 
symmetry protected topological phase characterized by 
the pair (G T pt, T Haldanc ) where G T pt is the symmetry 
group generated by time-reversal, parity, and translation 
symmetries. 

People have been using the degenerate spin-1/2 degrees 
of freedom at the ends of spin-1 chain, 51-53 and string or- 
der parameter 54 ' 55 to characterize Haldane phase. Both 
degenerate boundary spin-1/2 degrees of freedom and the 
string order parameter can be destroyed by perturba- 
tions that break spin rotation symmetry but have time- 
reversal, parity, and translation symmetries. So if wc 
used the degenerate spin-1/2 degrees of freedom and/or 
string order parameter to characterize the Haldanc phase, 
we would incorrectly conclude that the Haldane phase 
is unstable against those perturbations. Since the Hal- 
dane phase is indeed stable against those perturbations, 
we conclude that the (Gtpt ,T Haldanc ) characterization 
is more general. We also find that an arbitrary small 
perturbation that break parity symmetry will destabilize 
the Haldanc phase. This confirms a previously known 
result. 56 

In the first a few sections of this paper, we will review 
and introduce several tensor network renormalization ap- 
proaches. Then we will discuss new physical picture 




FIG. 1: A graphic representation of a tensor-network on a 
square lattice. A vertex represent a tensor 2ij/ e and the legs 
of a vertex carries the indices i, j, etc . Each link carries the 
same index. The indices on the internal links are summed 
over, which defines the tensor-trace. 

about phases and phase transitions emerged from the 
tensor network renormalization approach through several 
simple examples. 



II. A REVIEW OF A TENSOR NETWORK 
RENORMALIZATION METHOD 

In the following review, we will concentrate on 1+1D 
quantum systems. However, our discussion can be gener- 
alized to 2D statistical systems and the systems in higher 
dimensions. 

To describe the tensor renormalization approach intro- 
duced in Ref. 47, we first note that, after we discretize 
the space-time, the partition function represented by a 
space-time path integral Tre - ^ = J e~ J £ (in imagi- 
nary time) can be written as a tensor-trace over a tensor- 
network (see Fig. 1 for an example on 2D square lattice): 

Tre"^ = TjfeiThgjkTqkirTuts ■ ■ ■ = tTr ® 4 T. (2) 

ijkl--- 

where the indices of the tensor run from 1 to D (in other 
words the tensor T is a rank- four dimension- D tensor). 
Choosing different tensors T corresponds to choosing dif- 
ferent Lagrangian L. We see that calculating the tensor 
trace allows us to obtain properties of any physical sys- 
tems. 

Unfortunately, calculating the tensor-trace tTr is an 
NP hard problem in 1+1D and higher dimensions. 46 To 
solve this problem, Ref. 47 introduced an approximate 
real-space renormalization group approach which accel- 
erate the calculation exponentially. The basic idea is 
quite simple and is illustrated in Fig. 2. After finding 
the reduced tensor T", we can express tTr[T ® T ■ ■ ■ ] 
tTr[T" <g> T" ■ ■ ■ ] where the second tensor-trace only con- 
tains a quarter of the tensors in the first tensor-trace. 
We may repeat the procedure until there are only a few 
tensors in the tensor-trace. This allows us to reduce the 
exponential long calculation to a polynomial long calcu- 
lation. 

The actual implementation of the renormalization is a 
little more involved. 45 ' 47 For an uniform tensor-network 




(a) (b) (c) 



FIG. 2: The tensor T in the tensor network (a) has a dimen- 
sion D. After combine the two legs on each side into a single 
leg, the four linked tensors in (a) can be viewed as a single 
tensor T' in (b) with dimension D 2 . T' can be approximately 
reduced to a "smaller" tensor T" in (c) with dimension D and 
satisfies tTr[T' <g> T' ■ ■ ■] » tTr[T" ® T" ■ ■ ■]. 




FIG. 3: (Color online) (a) We represent the original rank-four 
tensor by two rank-three tensors, which is an approximate de- 
composition, (b) Summing over the indices around the square 
produces a single tensor T' . This step is exact. 

on 2D square lattice (see Fig 1), the renormalization 
group procedure can be realized in two steps. The first 
step is decomposing the rank-four tensor into two rank- 
three tensors. We do it in two different ways on the 
sublattice purple and green (see Fig. 4a). On purple sub- 

D 2 

lattice, we have T ru id = 2J s =i SiuisSadrs an d on green 

D 2 

sublattice, we have T ru id — zJs=i SzidsSirus- Note that 
r,l,u,d run over D values while s run over D 2 values. 
For such a range of s, the decomposition can always be 
exact. 

Next we try to reduce the range of s through an 
approximation. 47 Say, on purple sublattice, we view T ru id 




(a) (b) 1 'X ( C ) 



FIG. 4: (Color online) RG transformation of tensor-network 
produces a coarse grained tensor-network. 
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FIG. 5: (a) A CDL tensor, (b) A network of CDL tensor. 

as a matrix Mf u . rd = T TU u and do singular value decom- 
position 



M P = UAV^ A 



V 



A 2 
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(3) 



where Ai > A2 > A3 ... 
D cut singular values A s 



We then keep only the largest 



and define S; 



3 firs 



/atv; 



Si 



uls 



'XsUiu^. Thus, we can approximately express 



Truid by two rank-three tensors Si, S 3 



Truid 



D cut 



c 

Is^Sdrs- 



(4) 



8=1 



Similarly on green sublattice we may also define T ru id as 
a matrix Mf d . ru and do singular value decompositions, 
keep the largest D cut singular values and approximately 
express T ru id by two rank-three tensors S2, S4. 



T, 



D cut 

E 

7=1 



(5) 



After such decompositions, the square lattice is de- 
formed into the form in Fig. 4b. The second step is sim- 
ply contract the square and get a new tensor T' on the 
coarse grained lattice (see Fig. 3b). The dimension for 
the reduced tensor T" is only D cut which can be chosen 
to be D or some other values. This defines the (approxi- 
mated) renormalization flow of the tensor network based 
on the SVD. Let us call such an approach SVD based 
tensor renormalization group (SVDTRG) approach. 



III. THE PHYSICAL MEANING OF 
FIXED-POINT TENSORS 

Why the SVDTRG flow fails to produce a isolated 
fixed-point tensors? Levin has pointed out that one large 
class of tensors - corner double-line (CDL) tensors - are 
fixed points of the SVDTRG flow. 47 ' 57 A CDL tensor has 
a form (see Fig. 5a) 

T ruld (M\M 2 ,M 3 ,M 4 )=T ri 



r 2 ,u 2 ui l 2 l\ ,d\d 2 



= M 1 M 2 , 



(6) 



FIG. 6: The fixed-point tensor and its tensor network for 
topological phases. 



where the pair (ri,r2) represents the index r, etc . Such 
type of tensor (assuming \[T) = integer and n etc run 
from 1 to \TD) has exact decomposition 

-Dcut Dcut 

Truid — ^ S 3drs Sn us — StldsSirus 

8=1 S=l 

with I? cu t = D. This makes the CDL tensor to be the 
fixed-point tensor under the SVDTRG flow at least when 
M 1 = M 2 = M 3 = M 4 . If the initial tensor is changed 
slightly, the tensor may flow to a slight different CDL 
tensor, and hence the fixed point is not isolated. 

A tensor network of CDL tensor is plotted in Fig. 5b. 
If the system described by a tensor network does flow to 
a tensor network of CDL tensor, what does this means? 
From Fig. 5b, we see that, at the fixed point described 
by the CDL tensor, the degrees of freedom form clusters 
that correspond to the disconnected squares. Only the 
degrees of freedom within the same square interact and 
there are no interactions between different clusters. Such 
a fixed-point tensor (ie Lagrangian) describes a system 
with no interaction between neighbors and has a finite 
energy gap. So the system is in a trivial phase without 
symmetry breaking and without topological order (in the 
sense of long range entanglement). The li> J phase of 
the transverse-field Ising model 



-falsing = - ^(JviV?+i + hal) 



(7) 



is such a phase. We can check that the tensor network 
that describes the above transverse-field Ising model in- 
deed flows to a network of CDL tensor when fe> J. 

On the other hand, the low energy physics of such 
phases is really trivial. J — and h = 00 limit, the 
partition function of the transverse-field Ising model is 
described by a dimension-one tensor 



rpTRl 

1 ruld 



5r.l5u,\8l i8d,l 



(8) 



Thus we expect that the fixed-point tensor for the h^> J 
phase should be the above trivial dimension-one ten- 
sor. Unfortunately, the renormalization flow generated 
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by SVDTRG fails to produce such simple fixed-point ten- 
sor. 

The h <C J phase of the transverse-field Ising model is 
a symmetry breaking phase with two degenerate ground 
states. Under the SVDTRG flow, the tensor network 
of the transverse-field Ising model will flow to a tensor 
network described by a tensor of the form 

T = T cdl © T cdl (9) 

where T c di is a CDL tensor. On the other hand, the sym- 
metry breaking phase should be described by a simpler 
tensor eqn. (9) where T cc ii is the dimension-one tensor 
T TRI given by eqn. (8). 

However, the SVDTRG flow does a remarkable job of 
squeezing all the complicated entanglement into a sin- 
gle plaquette! Note that the initial tensor network may 
contain entanglement between lattice sites separated by 
a distance less than the correlation length. A CDL ten- 
sor contains only entanglement between nearest-neighbor 
lattice sites, while the ideal fixed-point tensor eqn. (8) has 
no entanglement at all (not even with next neighbors). 

What should be the fixed-point tensor for topological 
phases under SVDTRG flow? The essence of topological 
phases is pattern of long-range entanglement. 8-10 Since 
the CDL tensor has no long range entanglement, its can- 
not be used to describe topologically ordered phases. We 
believe that topological phases are described by fixed- 
point tensors of the following form (see Fig. 6) 

T = T clU ® T top (10) 

where T c di is CDL tensor and T top is a tensor that can- 
not be decomposed further as T top — T' cdl <g> T' t . Here 
T top will be the fixed-point tensor that characterizes the 
topological order. 



IV. FILTERING OUT THE LOCAL 
ENTANGLEMENTS 

We have seen that the SVDTRG flow can produce 
a network of CDL tensor which squeezes the non- 
topological short range entanglement into single plaque- 
tte. To obtain ideal fixed-point tensor (such as the one 
in eqn. (8)) we need to go one step further to remove 
the corner entanglement represented by the CDL tensor, 
ie we need to reduce corner matrix M z in the CDL tensor 
(see eqn. (8)) to a simpler form 

M i ab -> * 0il a 6il . (ii) 

In the following, we will describe one way to do so. 

In this approach, we modify the SVDTRG trans- 
formation described by Fig. 4 to a new transfor- 
mation described by Fig. 7. We will call the new 
approach tensor-entanglement-filtering renormalization 
(TEFR) approach. In the TEFR transformation, we first 
deform the square lattice Fig. 7(a) into the lattice Fig. 




(d) (c) 




(e) 



FIG. 7: (Color online) The TEFR procedure. We first de- 
form the tensor network on square lattice (a) to the tensor 
network on lattice (b) using the transformations in Fig. 3(a). 
We then replace the four rank-three tensors Si , S2 , S3 , Si in 

(b) by a new set of simpler rank-three tensors S[, S' 2 , S3, S4 
in (c) while keeping the tensor-trace of the tensor network 
approximately the same (see Fig. 8). The index on a solid 
link has a dimension D. The index on a dashed link has a 
dimension D' . Thus the tensor Si has a dimension D x D x D 
while the tensor S{ has a dimension D x D x D' . The SVD 
method is used to go from (a) to (b) [see Fig. 3(a)] and from 

(c) to (d) (see Fig. 9). The four S" on the small squares in 

(d) are traced over to obtain T' in (e). 



7(b) by expressing the rank-four tensor T in terms of two 
rank-three tensors Si, S3 or S2,Si [see Fig. 3(a)]. Then 
we replace the four rank-three tensors Si, S2, S3, 64 by 
a new set of simpler rank-three tensors S[ , S' 2 , S' 3 , S4 so 
that the tensor trace of S'i, S' 2 , S3, S4 (approximately) re- 
produces that of Si, S2, S3, S4 as shown in Fig. 8. In this 
step, we try to reduce the dimensions of the rank-three 
tensors, from D x D x D of Si to D x D x D' of S- with 
D' < D. If T is a CDL tensor, Si,S 2 ,S 3 ,S 4 will also 
be CDL tensors [see Fig. 10(a)]. In this case, the above 
procedure will result in a new set of tensors S[, S' 2 , S' 3 , S' 4 
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FIG. 8: (Color online) A tensor trace of a product of four 
rank-three tensors Si, S2, S3, £4 give rise to one rank- four ten- 
sor. Such a rank-four tensor can also be produced (approxi- 
mately) by a different set of rank-three tensors S[, S'2, S3, 1S4. 
We may choose S'i trying to minimize the dimension of the 
index on the dashed links. Or we may choose S[ trying to 
optimize some other quantities. 




(a) (b) 



FIG. 9: (Color online) The SVD method is used to implement 
the above two transformations. 



with a lower dimension as shown in Fig. 10(b). (For de- 
tails, see Appendix 1.) We then use the SVD (see Fig. 
9) to deform the lattice in Fig. 7(c) to that in Fig. 7(d). 
Last we trace out the indices on the small squares in Fig. 
7(d) to produce new tensor network on a coarse grained 
square lattice [see 7(e)] with a new tensor T'. The ten- 
sor T' will have a lower local entanglement than T. In 
fact, when T is a CDL tensor, the combined procedures 
described in Fig. 4 and Fig. 7 (see also Fig. 10 and 
11) can reduce T to a dimension-one tensor. Appendix 
1 contains a more detailed description of TEFR method. 



V. EXAMPLES 

In the section, we will apply the TEFR method to 
study several simple well known statistical and quantum 
systems to test the TEFR approach. We will calculate 
the free energy (for statistical systems) or the ground- 
state energy (for quantum systems). The singularities in 
free-energy /ground-state-energy in the large system limit 
will indicate the presence of phase transitions. We can 




(a) (b) 



FIG. 10: TEFR transformation of a CDL tensor. Here the 
index on a dash- line has a dimension one. 




(e) (d) 



FIG. 11: A tensor network of CDL tensor can be reduced to 
a tensor network of dimension-one tensor, (a) — > (b) uses the 
transformation in Fig. 10(a). (b) — » (c) uses the transforma- 
tion in Fig. 10(b). (c) — * (d) uses the transformation in Fig. 
9. (d) — » (e) uses the transformation in Fig. 4(b)— >(c). Here 
the index on a dash- line has a dimension one. 
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FIG. 12: (Color online) The temperature dependence of free 

energy F(T) and F" = f° r the 2D statistical Ising model. 
A phase transition can be seen around the temperature T « 
2.3. The vertical line marks the exact T c — 2.26919. 

also calculate the fixed-point tensors under the TEFR 
flow for each phase. The fixed-point tensors will allow 
us to identify those phases as symmetry breaking and/or 
topological phases. If we find a continuous phase transi- 
tion in the above calculation, we can also use the fixed- 
point tensor at the critical point to calculate the spec- 
trum of the scaling dimensions hi and the central charge 
c for the critical point (see Appendix 2). 



A. 2D statistical Ising model 

Our first example is the 2D statistical Ising model on 
square lattice described by 

H = -^<j i a j , <n = ±l (12) 
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FIG. 13: The graphic representations of (a) ^ T ruru , (b) 

2^2 ru ldTruluTldrdi and (c) y~l rTI ?rf T ru ldTldru- 



where Y^u j\ sums over nearest neighbors. The partition 
function is given by Z = Tr(e _/3ff ). Such a partition 
function can be expressed as a tensor-trace over a tensor- 
network on a (different) square lattice where the tensor 
T is given by 



T 



Ising 
1,2,1,2 



-4/3 



T, 



Ising 



2,1,2,1 



-4/3 



rrJsing _ 4/3 

J 1,1,1,1 — e ' 
others = 1. 



^Ising _ 4/3 
J 2.2,2,2 — e ' 



(13) 



Note that the spins are located on the links of the tensor- 
network (ie i in (12) labels the links of the tensor- 
network) . 

The statistical Ising model has a Z2 spin flip symme- 
try. Such a symmetry implies that the tensor is invariant 
under the following Z2 transformation 



T, 



Ising 
ruld 



= W rr ,W uu ,W lv W dd .Tlf*l dl , 



(14) 



where 



W 



1 

1 



We performed 14 steps of TEFR iteration described 
in Fig. 7 with D cut — 32. The calculated free energy 
F = -/3" 1 log(Z) = -Tlog(Z) is presented in Fig. 12. 
For temperature T > 2.38, the fixed-point tensor is a 
dimension-one tensor 



rpTRl _ 1 

J im — x > 



others = 0. 



(15) 



For temperature T < 2.19, the fixed-point tensor is a 
dimension-two tensor 



-1 

-'mi — L > 



1 

J 2222 — 1 ' 



others = 0. 



(16) 



For temperature 2.22 < T < 2.32, the fixed-point tensor 
is complicated with many non-zero elements. (We like to 
point out that the above fixed-point tensors are actually 
invariant fixed-point tensors. For details, see Appendix 

To visualize the structure of fixed-point tensors more 
quantitatively in different regions, we introduce the fol- 
lowing two quantities (see Fig. 13) 



Xi = 



C^2ru Truru) 
^^ruld ^TuluTldrd 



Xo = 



(Y2ru Tr 



Y2ruld TruldTldru 



(17) 
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FIG. 14: (Color online) The temperature dependence of Xi, 
X2, and the central charge c for the 2D statistical Ising model. 
The vertical line marks the exact T c = 2.26919. 



for the fixed-point tensors. Note that X\ and X2 is inde- 
pendent of the scale of the tensor T — * TT. We plotted 
X\, X2 and the central charge obtained from (53) for the 
fixed-point tensors Tj nv at different temperatures in Fig. 
14. (Note that r = t\ = i here.) We see that the differ- 
ent values of X\ and X2 mark the different phases, and 
the steps in X\ and X2, as well as the peak in c mark 
the point of continuous transition. 

From the structure of the fixed-point tensor, we see 
that the low-temperature phase T < 2.19 is a Z2 sym- 
metry breaking phase (where the fixed-point tensor has 
a form T Z2 = T™ © T TRI ) and the high- temperature 
phase T > 2.38 is a trivial disordered phase. Since 
dF/dT appear to be continuous and d 2 F/dT 2 appears 
to have an algebraic divergence at the transition point, 
all those suggest that the phase transition is a continuous 
transition, which is a totally expected result. From the 
Onsager's solution 58 sinh(2/T c ) = 1, one obtains the ex- 
act critical temperature to be T c — — — — r = 2.26919, 

. . . lo g(^2 +1 ) 

which is consistent with our numerical result. 

There is another interpretation of the fixed-point ten- 
sors. The fixed-point in the low temperature phase 
T z-2 = t tri t tri happen to be the initial tensor (13) 
in the zero temperature limit (3 — 00. Thus the flow of 
the tensor in the low temperature phase can be viewed 
as a flow of the temperature [3 — > 00. Using the same 
reasoning, we expect that the flow of the tensor in the 
high temperature phase can also be viewed as a flow of 
temperature (3 — > in the opposite direction. This would 
suggest that the fixed-point tensor should have a form 



T, 



TRI 

ild 



r, u,l,d — 1,2 



(18) 



which is very different from the trivial dimension-one ten- 
sor T™ in (15). We will explain below that the two 
tensors T TRI and J 1 ™ are equivalent. 

To see the equivalence, we would like to point out that 
two tensors T and T' related by 



tL 



(A )m(B )u'uT r ui d J 



(19) 



9 



give rise to same tensor trace for any square-lattice ten- 
sor network. In particular, if the singular values in the 
SVD decomposition in Fig. 3(a) has degeneracies, the 
coarse grained tensor will have the ambiguity described 
by the above transformation with A and B being orthog- 
onal matrices. Physically, the transformation (19) corre- 
sponds to a field redefinition. One can check explicitly 
that 



pTRI 

- ruld 



R r r> R U u' T^'un'd' Rl'lRd'd 



where 



R = -( 1 1 



1 

V2 V 1 



(20) 



Thus, T TRI and T TRI are equivalent under a field redef- 
inition. 

We have seen that the coarse grained tensors after the 
TEFR transformations may contain a field redefinition 
ambiguity (19). As a result, the Z 2 transformation on 
the coarse grained lattice may be different from that on 
the original lattice. The Z 2 transformation on the coarse 
grained lattice is related to the original Z 2 transformation 
by a field redefinition transformation. 

However, one can choose a basis on the coarse grained 
lattice by doing a proper field-redefinition transformation 
such that the Z 2 spin flip transformation has the same 
form on the coarse grained lattice as that on the original 
lattice. In such a basis, the disordered and the Z 2 sym- 
metry breaking phases are described by (G sym ,Ti nv ) = 
(Z 2 ,f™) and (G sym ,T inv ) = (Z 2 ,T Z »), respectively, 
where the Z 2 transformation is given by (14). This ex- 
ample shows how to use the pair (G sym , Ti nv ) to describe 
the two phases of the statistical Ising model. 

Next we like to use the fixed-point tensor Tj nv at the 
critical point to calculate the critical properties of the 
above continuous transition. We can use (49) to obtain 

(i) 

the invariant fixed-point tensor TV even off the critical 

point. From we calculate the matrices M ud , M lr , 
M ldru , and M lurd from (50), and obtain the eigenvalues 
\ n of those matrices. Using those eigenvalues, we can 
calculate the central charge and the scaling dimensions 
from (53), where r = n = i for the Ising model discussed 
here. By choosing Z? cu t = 64 and doing 9 iterations of 
Fig. 7 (that correspond to a system with 1024 spins) at 



the critical temperature T c 



log (72+1 



■, we find that 



c 

0.49942 
1/2 



hi 

0.12504 
1/8 



h 2 

0.99996 
1 



h 3 

1.12256 

9/8 



hi 

1.12403 

9/8 



where we also listed the exact values of the central charge 
and the scaling dimensions. We would like to mention 
that the TEFR approach does not suppose to work well 
at the critical point where the truncation error of keep- 
ing only a finite D cut singular values is large. It turns 
out that for Z? cu t = 64, the estimated relative truncation 



error is less than 10 -5 through 9 iterations. This allows 
us to obtain quite accurate central charge and scaling 
dimensions. The truncation error will grow for more it- 
erations. As a result, the agreement with the exact result 
will worsen. 

From this example, we see that the TEFR approach is 
an effective and efficient way to study phases and phase 
transitions. The TEFR approach allow us to identify 
spontaneous symmetry breaking from the structure of 
fixed-point tensors. The critical properties for the con- 
tinuous phase transition can also be obtained from the 
fixed-point tensor at the critical point. 



B. 2D statistical loop-gas model 

Next we consider a 2D statistical loop-gas model on 
square lattice. The loop-gas model can be viewed as a 
Ising model with spins on the links of the square lattice. 
However, the allowed spin configurations must satisfies 
the following hard constraint: the number of up-spins 
next each vertex must be even. The energy of spin con- 
figuration is given by 



(21) 



The partition function Z — Tr(e~@ H ) can be expressed as 
a tensor-trace over a tensor-network on the square lattice 
where the tensor T is given by 



T}2d = e /3(2d_3)/2 e /3(2tI_3)/2 e /3(2r_3)/2 e /3(2 ^ 3)/2 x 



$d+u+r+l=c 



(22) 



where the indices take the values r, u,l,d= 1,2 and 

8d+u+r+l=even = 1 if d + u + r + 1 — even, others = 0. 

To use the TEFR approach to calculate the partition 
function, it is very important to implement the TEFR in 
such a way that the closed-loop condition 

Truid — 0, when r + u + I + d = odd 
Srui — 0, when r + u + I = odd (23) 

are satisfied at every step of iteration. 

We performed 14 steps of TEFR iteration described 
in Fig. 7 with D cut = 32. Fig. 15 shows the resulting 
Xi, X 2 , and the central charge c = -log(Ao) obtained 
from (53) for the invariant fixed-point tensors (Ti nv ) ru id 
at different temperatures. We find that the low tem- 
perature phase is described by the trivial dimension-one 
fixed-point tensor T TRI (8) as indicated by X\ = X 2 = 1 . 
The high temperature phase has a non-trivial fixed-point 
tensor since X\ = X 2 = 2. We find that the invariant 
fixed-point tensor to have a form 



1 ruld 



if I + r + u + d = even, r, u,l,d= 1,2 



others = 0. 



(24) 
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FIG. 15: (Color online) The temperature dependence of Xi, 
X2, and the central charge c for the 2D statistical loop gas 
model. The vertical line marks the exact T c = 2.26919. 



The two fixed-point tensors T TRI and T LL correspond 
to the zero-temperature and infinite-temperature limit of 
the initial tensor T LG in (22). Both fixed-point tensors 
satisfy the closed- loop condition (23). 

It is well known that the high temperature phase of the 
statistical loop gas model is dual to the low temperature 
phase of the statistical Ising model, and the low temper- 
ature phase of the statistical loop gas model is dual to 
the high temperature phase of the statistical Ising model. 
The fixed-point tensors of the two models are related di- 
rectly 

rpLL p— 1 p— 1 p p 

1 ruld — ri rr' rl uu' 1 r>u'l'd' ri l'l n/ <l'd 

rpTKl _ d-1 d-1 'JtTRI p p 
1 ruld — n rr' tt uu' J -r'u'Vd' n -Vl tl d'd 

where R is given in (20). The critical properties at the 
continuous transition are also identical. 

One may wonder, if T LL and T Z2 are related by a field 
redefinition, then why the TEFR flow produces T LL as 
the fixed-point tensor of the high temperature phase of 
the loop gas, rather than T Z2 1 Indeed, if we implement 
the TEFR flow in the space of generic tensors, then both 
T LL and T Zq - , together with many other tensors related 
by some field redefinitions will appear as the fixed-point 
tensors of the high temperature phase of the loop gas. 
However, if we implement the TEFR flow in the space of 
tensors that the describe closed loops (ie satisfying (23)), 
then only T LL will appear as the fixed-point tensor. 



C. Generalized 2D statistical loop-gas model 

Third, let us discuss a more general 2D statistical loop- 
gas model described by a tensor network with tensor 



r GLG 

ruld 



-(r+l+u+d-6)f3+V5 rlud r 
e u d+u+r+l— e 

where r,u,l,d= 1,2 and 5 r i u d is given by 

£1111 = 62222 = 1, others = 0. 



(25) 



FIG. 16: (Color online) The V dependence of logz = 
log(Z)/iV and (logz)" = d 2 logz/dV 2 for the generalized 2D 
loop gas model. The sharpest peak for (logz)" and the logz 
currve are for a tensor network with TV = 4 7 tensors. Other 
peaks are for tensor networks with N — 4 6 , N — 4 5 , N = 4 4 , 
and N = 4 3 tensors. 



We note that when V = 0, the generalized loop-gas model 
become the one discussed above: T LG = T GLG \v=o- 

Let us consider the phases of the generalized loop-gas 
model along the (3 = line. In this case, the loop-gas 
model also has a Z2 symmetry: the tensor T GLG is in- 
variant under switching the 1 and 2 index: 



T 



GLG 



ruld 



= W r 



'W uu ,W w W dd ,T G ^, d ,, 



W x 



1 

1 



(26) 



We will call such a symmetry Zf symmetry. The closed- 
loop condition (23) can also be represented as a Z2 sym- 
metry of the tensor: 



t^2 = w rr ,w u 



'W u ,W ddl T r G ^, d ,, 



w z = 



1 
-1 



(27) 



The second Z2 symmetry will be called ivf symmetry. So 
along the (3 = line, the model has a x Zf symmetry. 

We performed the TEFR calculation for the general- 
ized loop-gas model along the (3 = line with Z3 cu t = 32. 
Fig. 16 shows the V dependence of the resulting log 
partition function per tensor, logz = log(Z)/N , and 
d 2 logz/dV 2 . Fig. 17 and Fig. 18 show the V depen- 
dence of Xi, X2, and central charge c for the resulting 
fixed-point tensor. We see that there is a phase transition 
at V c = 1.0985. The fixed-point tensor for the V < V c 
phase is found to be T LL and the fixed-point tensor for 
the V > V c phase is found to be T Z2 . The central charge 
is zero for the two phases on the two sides of the transi- 
tion. Thus the both phases have short range correlations. 
Using the (G sym , Ti nv ) notation, the two phases are char- 
acterized by (Zf x Z|,T LL ) and (Zf x Z|,T Za ). 

Since dlog(Z)/dV has no discontinuous jump at V c , 
the transition is a continuous phase transition. Also the 
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FIG. 17: (Color online) The V dependence of Xi, X 2 for the 
generalized 2D loop gas model. 
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FIG. 18: (Color online) The V dependence of the central 
charge c for the generalized 2D loop gas model. The sharpest 
peak is for a tensor network with N = 4 7 tensors. Other 
peaks are for tensor networks with N = 4 6 , N = 4 5 , N = 4 4 , 
and iV = 4 3 tensors. 

central charge is non zero at the transition point which 
suggests that the correlation length diverges at V c . This 
also implies that the transition is a continuous phase 
transition. From Fig. 18, we see that the transition is a 
central charge c = 1 critical point. 

We note that the two fixed-point tensors T LL and T 7,2 
give rise to the same X\ and X^. In fact, T LL and T Z2 
are related by a field-redefinition transformation (19): 

rpLL O— 1 E>— 1 rpZ2 T> T> 

1 ruld — n rr' n uu' 1 r'u'Vd' lt l'l lt <l'd 

where R is given in (20). Since X\ and Xi are invariant 
under the above transformation, so they are the same in 
the two phases. 

The generalized loop-gas model provides us an inter- 
esting example that the fixed-point tensors T LL and T 7,2 
that describe the two phases are related by a field re- 
definition transformation. In this case, one may won- 
der should we view the V < V c and V > V c phases as 
the same phase? A simple direct answer to the above 
question is no. It is incorrect to just use a fixed-point 



tensors to characterize a phase. We should use the 
pair (G sym ,Ti nv ) to characterize a phase. If we use the 
(G sym , Ti nv ) notation, the two phases are more correctly 
described by (Zf x Zf,T LL ) and (Zf x Zf,T Z2 ). The 
above field-redefinition transformation transforms T LL to 
T 7,2 . It also transform the symmetry transformations 
Zf x Zf to a different form Zf x Zf. Thus the pair 
(Zf x Zf,T LL ) is transformed to (Zf x Zf,T Z2 ) un- 
der the above field-redefinition transformation. There- 
fore, we cannot transform the pair (Zf x Zf,T LL ) to 
(Zf x Zf,T z ) using any field redefinition transforma- 
tion. This implies that the two phases (Zf x Zf,T LL ) 
and (Zf x Zf,T Z2 ) are different if we do not break the 
Zf x Zf symmetry. If we do break the Zf x Zf sym- 
metry, then the two phases are characterized by (1, T LL ) 
and (1,T Z2 ), where 1 represents the trivial group with 
only identity implying that the system has no symme- 
try. Without any symmetry, the two labels (1,T LL ) and 
{l,T z ' 2 ) are related under a field redefinition transfor- 
mation. This implies that the two phases (1,T LL ) and 
(l,T Z2 ) are the same phase if the system has no sym- 
metry. In Appendix 3, we will give a more detailed and 
general discuss about the definition of phases and its re- 
lation to the symmetry of the system. 

D. Quantum spin-1/2 chain and phase transition 
beyond symmetry breaking paradigm 

The generalized loop-gas model is closely related to the 
following ID quantum spin-1/2 model: 

H = Yy-^at +1 -Jatat +1 -hat\ (28) 

i 

Note that — of )/2 = J2i n i m °d 2 is a conserved 
quantity and the above model can be viewed as a ID 
hard-core boson model with next neighbor interaction. 
Such a ID quantum model can be simulated by a 2D 
statistical loop-gas model described by a tensor network 
with tensor T GLG (25). The loops in the loop gas cor- 
respond to the space-time trajectories of the hard-core 
bosons. The j3 = limit of the loop-gas model corre- 
spond to h = limit of the ID quantum model, and V 
corresponds to J. 

In the following, we will concentrate on the h = 
case. When h = 0, the Hamiltonian (28) has a Z2 x 
Zf symmetry, where Zf is the Z2 group generated by 
cr x and Zf is the Z2 group generated by <r z . Such a 
Zf x Zf symmetry corresponds to the Zf x Zf symmetry 
discussed in last section. 

When J < J c = 1 , the system is in a symmetry break- 
ing phase. We will use a pair of groups (G sym , G grn d) 
to label such a symmetry breaking phase, where G sym 
is the symmetry group of the Hamiltonian and G grn d is 
the symmetry group of a ground state. So the J < J c 
phase of our system is labeled by (G sym , G grn d) = (Zf x 
Zf , Zf ), where the ground state breaks the Zf symmetry 
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(a) The graphic representation of 
(b) e~ STHi can be viewed as a 



but has the Z% symmetry. When J > J c , our system is 
in a different symmetry breaking phase, which is labeled 
by (G sym , Ggmd) = (Z% x Zf,Z|). 

Under the correspondence between the loop gas and 
our quantum spin- 1/2 model, the V < V c and V > V c 
phases of the loop-gas model discussed above will corre- 
spond to the above two Z 2 symmetry breaking phases for 
J < J c = 1 and J > J c respectively. The results from 
the loop gas imply that the two Z 2 symmetry breaking 
phases are connected by a continuous phase transition. 
Such a continuous phase transition can be viewed as two 
Z 2 symmetry breaking transition happening at the same 
point. 

Strictly speaking, the above continuous transition does 
not fit within the standard symmetry breaking paradigm 
for continuous phase transitions. According to Landau 
symmetry breaking theory, a continuous phase transi- 
tion can only happen between two phases labeled by 
(G sym ,G grnd ) and (G sym , G' glmd ) where G grnd is a sub- 
group of G grnd or Gg rnd is a subgroup of G grn d- We see 
that the continuous transition between the (Z 2 x Z 2 , Z 2 ) 
phase and (Z 2 x Z 2 ,Z 2 ) phase does not fit within the 
standard Landau symmetry breaking theory. Such a 
continuous transition is another example of continuous 
transitions that is beyond the Landau symmetry break- 
ing paradigm. 29 ' 59-6 Even though the above continuous 
phase transition does not fit within the standard symme- 
try breaking theory, such a transition can still be studied 
through the tensor network renormalization approach. In 
particular, the tensor network renormalization approach 
reveals that the transition is a central charge c = 1 criti- 
cal point. 



E. Quantum spin-1 chain 

Last, we consider a quantum spin-1 chain at zero tem- 
perature. The Hamiltonian is given by 



H = ^2(Si- S l+1 + U{S*f + BS* 



(29) 



The model has a translation symmetry, a Z 2 time rever- 
sal symmetry S y — ► — S y and Z 2 parity symmetry. The 
imaginary-time path-integral of the model can be written 



Tre-™ = lim Tr( e -«Tff. e -tfrH.jM 

M—+oo 



where St = t/M and 

H e = (Si-S i+1 + U(S!) 2 + ^(Sf + Sf +1 )) 

i—even 

H =Y J (Si-S i+1 + U(S?) 2 + ^(Sf + Sf +1 ) 



i— odd 



Note that H e and H Q are sum of non overlapping terms 
and thus 



-StH, 



•= n 

<,=evei 

n 



c -8rH 0= 1{ , 
i— odd 



-SrHi 



Hi in the above expression has a form 



II - S.-8 H ('(S?) 2 + Z(S! + Sf +1 ) 



The matrix elements of e 5rHi can be represented by a 
rank-four tensor (see Fig. 19) 

where mj,m^ label of the quantum states of SV There- 
fore the imaginary-time path-integral Tre _Tff can be ex- 
pressed as a tensor trace over a tensor network of T (see 
Fig. 20(a)) 

Tre" rH = tTrOi T. 

Now we can use the TEFR approach to evaluate the ten- 
sor trace and to obtain the fixed-point tensors. This al- 
lows us to find the zero-temperature phase diagram and 
the quantum phase transitions of the spin-1 model (29). 

However, for small St, the tensor network is very 
anisotropic with very different correlation lengths (mea- 
sured by the distance of the tensor network) in the time 
and the space directions. This causes large truncation er- 
rors in the TEFR steps Fig. 7(a)-»(b) and Fig. 7(c)->(d) 
[or Fig. 3(a)]. We need to perform anisotropic coarse 
graining to make the system more or less isotropic be- 
fore performing the TEFR coarse graining. One such 
anisotropic coarse graining in time direction is described 
in Fig. 20. Each coarse graining step of Fig. 20 triple 
the effective St. After making the effective St to be of 
order 1, we switch to the isotropic coarse graining in Fig. 
7. 

We performed 11 steps of TEFR iteration in Fig. 2 
(which is implemented as 22 steps of TEFR iteration in 
Fig. 7) with D cut = 37 for various value of U and B. Af- 
ter the TEFR iteration, a initial tensor T(U, B) flows to 
a fixed-point tensor. As we vary U and B, the different 
initial tensors may flow to different fixed-point tensors 
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FIG. 20: (Color online) A T a -T b tensor network of size L x M 
is reduced to a T^-T(, tensor network of size L x M/2. From 
(a) to (b), the local deformation Fig. 3(a) is used. From (b) 
to (c), the local deformation Fig. 9 is used. From (c) to (d), 
a local deformation similar to that in Fig. 3(b) is used. 



which represent different phases. A more detailed de- 
scription of the renormalization flow of the tensors can 
be found in Appendix 4. 

To quantitatively plot the fixed-point tensors, we can 
use X\, X2 (see (17)) , and A SU m- Here A sum is given by 
Asn 



matrix M? urd 



A s /Ai| and Aj are the singular values of the 
obtained from the fixed-point tensor (see 
(50) and (3)). In Fig. 21, we plot the X\, X 2 , and A sum 
of the resulting fixed-point tensors for different value of 
U and B. 

We see that our spin-1 system has four different phases: 
TRI, Zf , Z y , and Haldane. The transition between the 
Haldane phase and the Z 2 phase happens at B c = 0.405 
along the U = line, which agrees with the Haldane gap 
A = 0.4097 = B c obtained in Ref. 65. The transition 
between the Haldane phase and the TRI phase happens 
at U c = 1.0 along the B = line, which agrees with the 
result U c = 0.99 obtained in Ref. 66. 

The fixed-point tensor in the TRI phase is the triv- 
ial dimension-one tensor T TRI with (Xi, X2, A sum ) = 
(1,1,1). So the TRI phase has no symmetry breaking. 
The fixed-point tensor in the Z\ and Z% phase, up to 
a field redefinition transformation, are the dimension- 
two tensor T Z2 which is a direct sum of two dimension- 
one tensor T™ with (Xi,X 2 , A sum ) = (2,2,2). So 
the Z\ and the iff phases are Z 2 symmetry breaking 
phases. The fixed-point tensor in the Haldane phase 
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FIG. 21: (Color online) The color-map plot of (r,g,b) = 
(Xi/3, X2/3, l/2Asum) shows the phase diagram of our spin-1 
system (29) which contains four phases, a trivial spin polar- 
ized phase "TRI", two Z2 symmetry breaking phases ".Zf" 
and "Zf", and a Haldane phase "Haldane". The transition 
between those phases are all continuous phase transitions. 
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FIG. 22: (Color online) The color-map plot of the central 
charge as a function of U, B. The color is given by (r,g,b) = 
(eg, eg, cio), where Ci is the central charge obtained at the i th 
iteration of Fig. 2. The lines with non-zero central charge 
are the critical lines with gapless excitations. The model (29) 
has a U(l) spin rotation symmetry along the U = line. The 
U — line between the two Z2 phases is a gapless state. 



is a dimension-four tensor. However, many different 
dimension-four tensors can appear as the fixed-point ten- 
sors. All those dimension-four fixed-point tensors have 
(Xi,X2, A sum ) = (1, 1,4). In fact they are all equivalent 
to the following CDL tensor 



T 



Haldane 



T(a v ,a v ,a y ,a y ) 



(31) 



through the field-redefinition transformation (19), where 
T(cr y ,a v ,a v ,(T y ) is given in (6). In fact, one can show 
that the following four CDL tensors, T(a x ,a x ,cr x , a x ), 
T(cr y , a y ,a y ,cr y ), T(a z 7 <j z ,cr z , a z ), and T(a°, er°, a , <r°), 
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are all equivalent through the field-redefinition transfor- 
mation. (Here a is the two by two identity matrix.) So 
they can all be viewed as the fixed-point tensors in the 
Haldane phase. Since the fixed-point tensor in the Hal- 
dane phase is a CDL tensor, the Haldane phase, like the 
TRI phase, also does not break any symmetry. 

From the fixed-point tensors, we can also obtain the 
low energy spectrum of our spin-1 system (see Appendix 
2). Using the fixed-point tensors of the four phases, we 
find that all four phases have a finite energy gap. We 
also find that the ground states of the TRI and Haldane 
phases are non degenerate while the ground states of the 
Zf and Zf phases have a two-fold degeneracy. This is 
consistent with our result that the TRI phase and the 
Haldane phase do not break any symmetry while the Zf 
and Zf phases break a Z2 symmetry. 

The central charge calculated from the resulting fixed- 
point tensors is plotted in Fig. 22 (see (53)). The lines 
of non-zero central charge in Fig. 22 mark the region of 
diverging correlation length, which correspond to line of 
continuous phase transition. The energy gap closes along 
those lines. Compare Fig. 22 with Fig. 21, we see that all 
the phase transitions in the phase diagram Fig. 21 cor- 
respond to those critical lines. Those phase transitions 
are all continuous phase transitions. The transitions be- 
tween Haldane and Zf , Haldane and Zf, TRI and Zf, 
are described by central charge c — 1/2 critical point. 
The transitions between Zf and Zf , Haldane and TRI, 
are described by central charge c = 1 critical point. 

The continuous phase transition between the TRI and 
Zi phases, and between the Haldane and Z2 phases are 
Z 2 symmetry breaking transitions. The phase transition 
between the Zf and the Zf phases is of the same type as 
that discussed in section VD along the h — line. Along 
the U = line, our spin-1 model has a 17(1) spin rotation 
symmetry around the S x axis. The gapless phase at B > 
B c and U = is close to a ?7(1) symmetry breaking phase 
where spins are mainly in the y-z plane with a finite and 
uniform S x component. When U ^ 0, the U(l) spin 
rotation symmetry is broken. When U > 0, the spin-1 
model is in the Zf phase with Sf cx (— )\ When U < 0, 
the spin-1 model is in the Zf phase with Sf oc (— )\ 

Both the TRI phase and the Haldane phase have the 
same symmetry, but yet they are separated by phase 
transitions. This suggests that the TRI phase and the 
Haldane phase are distinct phases. But we cannot use 
symmetry breaking to distinguish the two phases. So let 
us discuss Haldane phase in more detail. 



F. Haldane phase — a symmetry protected 
topological phase 

The first question is that weather the TRI phase and 
the Haldane phase are really different? (For a detail dis- 
cuss on the definitions of phases and phases transitions, 
see Appendix 3.) Can we find a way to deform the two 
phases into each without any phase transition? At first 



sight, we note that the fixed-point tensor for the Haldane 
phase are CDL tensors. According to the discussion in 
section III, this seems to suggest that the Haldane phase 
is the trivial TRI phase described by T TRI . 

In fact, the answer to the above question depend on 
the symmetry of the Hamiltonian. If we allow to deform 
the spin-1 Hamiltonian arbitrarily, then the TRI phase 
and the Haldane phase will belong to the same phase, 
since the two phases can be deformed into each other 
without phase transition. If we require the Hamiltonian 
to have certain symmetries, then the TRI phase and the 
Haldane phase are different phases (in the sense defined 
in Appendix 3). 

To demonstrate the above result more concretely, let us 
consider a spin-1 system with time-reversal, parity (spa- 
tial reflection), and translation symmetries. The time- 
reversal symmetry requires the Hamiltonian H and the 
corresponding tensor T (see 30) to be real. Since the real 
Hamiltonian must be symmetric, time-reversal symmet- 
ric tensor must also satisfy 

T ru ld = TdX ur (32) 

(Note that the time-reversal transformation exchanges 
(mj,mj + i) <-> (m-,m- +1 ) in (30)) The parity symmetry 
requires the tensor T to satisfy 

Truld — T urc i[. (33) 

(Note that the parity exchanges (mj,mQ <-» 
(m i+ i , m' i+1 ) in (30)) The tensor network considered 
here is always uniform which ensure the translation sym- 
metry. Our spin-1 Hamiltonian (29) has time-reversal, 
parity, and translation symmetries. 

Using our numerical TEFR calculation, we have 
checked that the Haldane phase is stable against any 
perturbations that respect the time-reversal, parity, and 
translation symmetries. In other words, we start with a 
tensor T that respects the those symmetries and flows to 
the fixed-point tensor T Halda,ne . We then add an arbitrary 
perturbation T — > T + ST that also respect the time- 
reversal, parity, and translation symmetries. We find 
that T+ST also flows to the fixed-point tensor T Haldane as 
long as ST is not too large. This result suggests that the 
time-reversal, parity, and translation symmetries protect 
the stability of the Haldane phase. This implies that in 
the presence of those symmetries, the TRI phase and the 
Haldane phase are always different. We cannot smoothly 
deform the Haldane phase into the TRI phase through 
the Hamiltonian that have the time-reversal, parity, and 
translation symmetries. 

If we add a perturbation ST that has time-reversal 
symmetry but not parity symmetry, our numerical re- 
sult shows that the perturbed tensor T + ST will fail to 
flow to T Haldanc . Instead, it will flow to the trivial fixed- 
point tensor T TRI . This implies that, without the parity 
symmetry, the Haldane phase is unstable and is the same 
as the trivial phase TRI, which agrees the result obtained 
previously in Ref. 56. 




To summarize, for Hamiltonian with time-reversal and 
parity symmetries, the Haldane phase and the TRI phase 
are different phases, despite both phases do not break 
any symmetries. So the difference between the Haldane 
phase and the TRI phase cannot be describe by Landau 
symmetry breaking theory. In this sense the Haldane 
phase behaves like the topologically ordered phases. 

On the other hand, for Hamiltonian without parity 
symmetry, the Haldane phase and the TRI phase are the 
same phase. We see that the Haldane phase exists as 
a distinct phase only when the Hamiltonian have time- 
reversal, parity, and translation symmetries. This behav- 
ior is very different from the topological phases which are 
stable against any perturbations, including those that 
break the parity symmetry. Therefore, we should re- 
fer Haldane phase as a symmetry protected topologi- 
cal phase. We like to stress that the Haldane phase is 
only one of many possible symmetry protected topolog- 
ical phases. Many other symmetry protected topologi- 
cal phases are studied using projected symmetry group 
in Ref. 17,29. For example, without symmetry, there is 
only one kind of Zi spin liquid. In the presence of spin 
rotation symmetry and the symmetries of square lattice, 
there are hundreds different Z^ spin liquids. 29 All those 
i?2 spin liquids can be viewed as symmetry protected 
topological phases. 

Since the time-reversal, parity, and translation sym- 
metries play such an important role in the very existence 
of the Haldane phase, it is not proper to characterize the 
Haldane phase just using the fixed-point tensor T Ral ' ia - ne . 
We should use a pair (Gtpt, T Haldanc ) to characterize the 
Haldane phase. Here Gtpt is the symmetry group gen- 
erated the time reversal, parity, and translation transfor- 
mations. Using the more complete notation, we can say 
that (Gtpt, T Haldanc ) and (G T pt,T™) correspond to 
different phases, while (G T t, T Haldanc ) and (G T t,T tr1 ) 
correspond to the same phase. Here Gtt is the symme- 
try group generated only by the time-reversal and trans- 
lation transformations. In terms of the new notation, the 
four phases, TRI, Haldane, and the two Z2, are labeled 
by (Gtpt,T™), (Gtpt, T Haldanc ) ; (G TPT ,T^). and 
(Gtpt, T 2 " 1 ). For a more detailed discussion on how the 
symmetry transformations act on the fixed-point tensor, 
see Appendix 5 

In this paper, we propose to use the pair 
(GTPT,T Kaldanc ) to characterize the Haldane phase - 
a symmetry protected topological phase. Such a nota- 
tion allows us to show that the Haldane phase is sta- 
ble against any perturbations that respect time-reversal, 
parity, and translation symmetry. We feel that the new 
characterization capture the essence of Haldane phase. 
People have proposed several other ways to characterize 
the Haldane phase. One of them is to use the emergent 
spin-1/2 boundary spins at the two ends of a segment of 
the isotropic spin-1 chain in the Haldane phase. How- 
ever, such a characterization fails to capture the essence 
of Haldane phase. The Haldane phase still exists in the 
presence of an uniform magnetic field: B Sf , while the 




FIG. 23: (Color online) The color-map plot of (r,g,b) = 
(Xi,X2, 1/Asum) shows the phase diagram of the spin-1 sys- 
tem (34) which contains four phases, a trivial phase "TRI", 
two Z2 symmetry breaking phases and "Zf ", an< l a Hal- 

dane phase "Haldane". The transition between those phases 
are all continuous phase transitions. 
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FIG. 24: (Color online) The central charge as a function of 
U, B. The lines with non-zero central charge are the critical 
lines with gapless excitations. The model (34) has a U(l) 
spin rotation symmetry along the U = line. The U = line 
between the two Z2 phases is a gapless state. 



four-fold degeneracy from the emergent spin-1/2 bound- 
ary spins is lifted by such an uniform magnetic field. In 
the presence of an uniform magnetic field, we can no 
longer use the degeneracy of the boundary spin-1/2 spins 
to detect the Haldane phase. 

People also proposed to use the string order parameter 



lim 

l»— *oo 



Sf 



^0 



to characterize the Haldane phase. But the string order 
parameter vanishes in the presence of perturbations that 



16 



U 








- 









B 



0.5 




B 



FIG. 25: (Color online) The color-map plot of (r,g,b) — 
(Xi/3, X2/3, l/2A S um) shows the phase diagram of the spin-1 
system (35) which contains four phases, a trivial spin polar- 
ized phase "TRI", two Z2 symmetry breaking phases ".Zf" 
and "Zf", and a Haldane phase "Haldane". The transition 
between those phases are all continuous phase transitions. 



contain odd number of spin operators, such as 
8H = Y,M(S z l ) 2 S! +1 + (S* +1 ) 2 St] 

i 

i 

This is because the total number of S z = +1 and 
S* — — 1 sites mod 2 is not conserved in the presence 
of the above perturbation. However, the vanishing string 
order parameter does not imply the instability of Hal- 
dane phase. The above perturbation respects the time- 
reversal and parity symmetries. The Haldane phase is 
stable against such a perturbation. Therefore, the string 
order parameter also fails to capture the essence of Hal- 
dane phase. 

To illustrate the stability of Haldane phase, we also cal- 
culated the phase diagrams of following two spin-1 mod- 
els. The first model is given by 



H 



E 



(Si ■ Si+i + U(S?) 2 + B(S* + Sf)) (34) 



The phase diagram (the color map plot of X\, X%, and 
Asum) and the central charge are plotted in Figs. 23 and 
24. The second model is given by 



H = J2 (Si ■ S i+1 + U(S*f) 



B 
~2 



E 



(Sf(Si + x) + S i+ 



z\2 



2Sf 



(35) 



The phase diagram and the central charge are plotted in 
Figs. 25 and 26. 

We see that in both models, the Haldane phase and 
the TRI phases are separated by phase transitions, indi- 
cating the stability of the Haldane phase. (Note that the 



FIG. 26: (Color online) The color-map plot of the central 
charge as a function of U, B. The color is given by (r,g,b) = 
(eg, eg, cio), where Ci is the central charge obtained at the i th 
iteration of Fig. 2. The lines with non-zero central charge are 
the critical lines with gapless excitations. 
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FIG. 27: (Color online) The fixed-point wave function of the 
Haldane phase can be obtained from the imaginary-time evo- 
lution of the fixe-point tensor network of the Haldane phase 
T Haidano = T {a y , a v , a v , a y ) , which is formed by the corner 
matrix a y . (ie the boundary of the tensor network gives raise 
to the fixed-point wave function.) The dimension two of the 
corner matrix implies that each index of the corner matrix is 
associated with a pseudo spin-1/2. Thus each coarse grained 
site has two spin-1/2 pseudo spins. The corner matrix a y ties 
the two pseudo spins on the neighboring sites into a pseudo 
spin singlets. Thus the fixed-point wave function of the Hal- 
dane phase can be viewed as a dimer state formed by nearest 
neighbor pseudo spins. 



stability of the Haldane phase is defined by the impos- 
sibility to deform the Haldane phase to the TRI phase 
without a phase transition within a symmetry class of 
Hamiltonian.) Both the models have time-reversal and 
parity symmetries. In particular, the second model has 
only time-reversal and parity symmetries. This suggests 
that the Haldane phase is stable if the Hamiltonian has 
time-reversal and parity symmetries. 

Before ending this section, we like to point out that 
the form of the fixe-point tensor of the Haldane phase 
T Haidanc _ T (ai> , a y , u v , a*) implies the following fixed- 
point wave function for the Haldane phase. Since the 
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FIG. 28: (Color online) The paradigm of topological order. 



fixed-point tensor T Haldano = T{a y , a y , &» , a y ) has a di- 
mension four, the effective degrees of freedom in the 
fixed-point wave function is four states per site [ie per 
coarse grained site). The four states can be viewed as 
two spin- 1/2 pseudo spins. The fixed-point wave func- 
tion is then formed by the pseudo spin singlet dimers 
between nearest neighbor bonds (see Fig. 27). 



VI. SUMMARY 

In this paper, we introduce a TEFR approach for sta- 
tistical and quantum systems. We can use the TEFR 
approach to calculate phase diagram with both symme- 
try breaking and topological phases. We can also use the 
TEFR approach to calculate critical properties of con- 
tinuous transitions between symmetry breaking and/or 
topological phases. The TEFR approach suggests a very 
general characterization of phases based on the symme- 
try group G sym and the fixed-point tensor T- mv . The 
(G sym , Ti nv ) characterization can describe both symme- 
try breaking and topological phases. Using such a char- 
acterization and through the stability of the fixed-point 
tensor, we show that the Haldane phase for spin-1 chain 
is stable against any perturbations that respect time- 
reversal, parity, and translation symmetry. This suggests 
that the Haldane phase is a symmetry protected topolog- 
ical phase. 

The TEFR approach introduced here can be applied to 
any strongly correlated systems in any dimensions. It has 
a potential to be used as a foundation to study topolog- 
ical phases and symmetry protected topological phases 
in any dimensions. In this paper, we only demonstrate 
some basic applications of the TEFR approach through 
some simple examples. We hope those results will lay the 
foundation for using the TEFR approach to study more 
complicated strongly correlated systems in 1+2 and in 
1+3 dimensions. 

To gain a better understanding of the tensor net- 
work renormalization approach, let us compare the tensor 
network renormalization approach with the block-spin 



renormalization approach. 24 Let us consider a spin sys- 
tem whose partition function can be represented by a 
tensor network of a tensor T(g), where g is a coupling 
constant of the spin system. After performing a coarse 
graining of the lattice, the system is described by a coarse 
grained tensor network of tensor T". The transforma- 
tion T(g) — > T" represent tensor network renormalization 
flow. In the block-spin renormalization approach, we in- 
sist on the coarse grained system have the same form as 
described by a new coupling constant g'. So we try to find 
a g' such that T(g') best approximates T' . This way, we 
obtain a block-spin renormalization flow g — > g' . We see 
that tensor network renormalization flow is a transfor- 
mation in the generic tensor space, while the block-spin 
renormalization flow is a transformation is the subspace 
T(g) parameterized by the coupling constant g. We be- 
lieve that it is this generic flow in the whole tensor space 
allows the tensor network renormalization approach to 
describe topological phases. 

In this paper, we stress that the fixed-point tensor of 
the tensor network renormalization flow can be used as 
a quantitative description of topological order. This is 
just one aspect of topological order. In fact, theory of 
topological order has a very rich structure and covers a 
wide area (see Fig. 28). Turaev-Viro invariant 67 on 3D 
manifold can be viewed as a tensor trace of the fixed- 
point tensors for string-net models. 12 Therefore, the ten- 
sor network renormalization approach should apply to 
all time reversal and parity symmetric topological orders 
in 1+2 dimensions which are described by the string-net 
condensations. 11 ' 12 This type of topological order include 
emergent gauge theory, Chern-Simons theory, and quan- 
tum gravity in 1+2 dimensions. 12 ' 68,69 The mathematical 
foundation that classifies this type of fixed-point tensors 
is the tensor category theory. It is natural to expect that 
the topological order in 1+3 dimensions and the asso- 
ciated emergent gauge bosons, emergent fermions, and 
emergent gravitons 34 ' 70 can also be studied through the 
tensor network renormalization approach. The tensor 
network renormalization approach also provides a new 
way to calculate lattice gauge theories. We see that 
physics is well inter connected. It appears that the ten- 
sor network renormalization is a glue that connects all 
different parts of physics. 



VII. ACKNOWLEDGEMENTS 



We would like to thank Michael Levin, Frank Ver- 
straete, and Michael Freedman for many very helpful 
discussions. This research was supported by the Founda- 
tional Questions Institute (FQXi) and NSF Grant DMR- 
0706078. 



18 



a \/ d 
b / ^ c 



FIG. 29: A rank-four tensor expressed as tensor network. 
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FIG. 30: The entanglement-filtering procedure. 



1. A detailed discussion of the TEFR algorithm 

The key step in the TEFR calculation is the entan- 
glement filtering operation from Fig. 7(b) to Fig. 7(c). 
In this section we will discuss several ways to implement 
such an entanglement filtering calculation. 



a. The linear algorithm 



To find those new described in Fig. 
minimize the following cost function: 



Ct 



ijkl 



we may 



(36) 



^llia^2ij/3S3jk-fS4klS 



ijkl 



ia^^ijP^^jk^^iklS 



where a, (3, 6, 7 = 1, . . . D and i, j, k,l = 1, . . . D' 

This cost function can be minimized by solving a sets 
of least square problems iteratively. For example, we 
may first fixed S' 2 , S' 3 , S' 4 and solve a least square problem 
for S[, then we fixed S[,S 3 ,S / 4: to solve a least square 
problem for S 2 , etc. In practise, to reduce the trapping 
by local minima, we may start with D' = D — 1 and 
gradually decrease D' . 



b. The nonlinear algorithm 

Although the linear algorithm is an efficient way to do 
the entanglement-filtering calculation in Fig. 8, however, 
this algorithm could not avoid trapping by local minima 
in generic case. To solve this problem, here we introduce 
a non-linear algorithm which can avoid the trapping by 
local minima. 

In the first implementation of TEFR discussed above, 
we choose S^ to minimize the dimension of the one of the 
index of the rank-three tensors S 4 '. The second imple- 
mentation of TEFR is still described by Fig. 7. But now 
we choose S£ to minimize the entanglement entropy on 
the diagonal links in Fig. 7(d). 

Let us define the entanglement entropy on a link. We 



1,3 



may view the tensor network in Fig. 29, as matrix M,, 
where i = (a, b) and j — (c, d). We then perform an 
SVD of M which produces the singular values A^ that 
are ordered as Ai > A 2 > ... > 0. We call A^ the singular 
values on the horizontal link in Fig. 29. The entangle- 
ment entropy on the horizontal link in Fig. 29 is defined 



as 5" — - J2i Aj log(Aj), where Ai = A;/X]A 
non-linear algorithm, we choose S[, S' 2l S' 3 
minimize the following cost function 



and 5*4 



In the 

to 



Cnl — AS' sum + (37) 

II E Si lia S2ijfjS3jk~/S4kl5 " E SuiaSz'ijpSzj^SiMs | 
ijkl ijkl 

where S snm is the sum of entanglement entropy on the 
four diagonal links in Fig. 7(d). Here A should be as 
small as possible. 

We note that when we reduce the entanglement en- 
tropy on the diagonal links in Fig. 7(d), the singular 
values on those diagonal links decay faster. This allows 
us to drop indices associated with the very small singular 
values and hence reduce the dimension on those diagonal 
links. When we use Fig. 3(b) to transform Fig. 7(d) to 
Fig. 7(e), the resulting tensor T' will have a small di- 
mension. This approach can again reduce a CDL tensor 
T to a trivial dimension-one tensor. 

The non-linear algorithm is a very general approach, 
but the calculation is very expensive. In the paper, we 
will not use this algorithm. 



c. The scaling-SVD algorithm for CDL tensors 

In this section, we introduce the third algorithm, which 
is very simple and efficient. We will call such algorithm 
scaling-SVD algorithm, which can be used to find S[, S' 2 , 
S' 3 , and S'4 in Fig. 8 with smaller dimensions. 

In the scaling-SVD algorithm, we introduce four diag- 
onal weighting matrices: 



.12 



w} 2 5u>: 



w lk' = wf5 kk , 



w 23 , 
n' 



wt l b~ w \ 



on those links i, j, k, I in Fig. 30 and initialize them as 1. 
We first perform the scaling-SVD calculation on the pair 
Si and S'2 by introducing a matrix: 



and do the SVD decomposition M 12 = UAV^: 
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We then update 5i, 52 as: 



A} 2 U al;i /Jw] 



5 



We can update 52, 53 in the same way, then we can fur- 
ther update 53, 54 and 54, S\. After one loop update for 
those 5j's, we update all the weighting factors w 12 etc as 
w 12 = A 12 etc on corresponding links. 

After several loop iterations, we truncate the smallest 
singular values of Aj 2 , A 23 , A 3 . 4 , A 41 and obtain new rank- 
three tensors S[, 5 2 , 5 3 , 5 4 with a smaller dimension D' < 
D on the inner links. If 5i, 5 2 , 53, and 54 are CDL 
tensor, we find: 

Si' Ha S2iji3S3j k ^S4 kls CX Sn ia S2ijfjS3j kl S4 k is 

after truncation. By proper normalizing those 5-, we can 
find the solution of 

^llia^ij0^3jk^^iklS ~ ^ ^Uia^ij/j^jk-y^iklS ~ 0- 

(38) 

If 5i, 52 , 53, and 54 are not CDL tensor, the truncated 
S[, 5 2 , 5 3 , and 5^ cannot satisfy (38) even after the 
rescaling of S[, 5 2 , 5 3 , and 5 4 . In this case, we reject the 
truncated S[, S' 2 , S' 3 , and 5 4 and the scaling-SVD fails 
to produce simplified tensors. To summarize, the scaling- 
SVD algorithm can simplify CDL tensors, but not other 
tensors. 

To understand why the above scaling-SVD algorithm 
can work for CDL tensors, we can consider the simplest 
CDL tensor T ru i d (M 1 , M 2 , M 3 , M 4 ) with corner matri- 
ces 



M 1 = M 2 = M 3 = M 4 



1 

x 



where < x < 1. 

See in Fig. 31, after one loop iteration, the 4 x 4 di- 
agonal weighting matrix w = u>i n ® wout on each link is 
updated as: 



wi % 



1 

x 



WQut 



1 

x 



Where In represents the inner line and Out represent 
the outer line of the square as show in Fig. 31. However, 
after n loop iteration, we find: 



Wlr, 



1 

x" 



WQut 



1 

x 



The striking thing is that the weighting matrices on the 
outer lines arc never scaled but the weighting matri- 
ces (as well as the singular values on the corresponding 
link) on the inner lines are scaled to its n power! After 
enough times of loop iterations, x n can be smaller than 
the machine error. The inner line can be removed and 
we successfully filter out those local entanglements in the 
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FIG. 31: For the simplest CDL tensor, the outer line singular 
values are not scaled after one loop iteration, however, the 
inner line singular values are scaled to a power. This nontriv- 
ial scaling behavior is the key why the scaling-SVD algorithm 
can filter out those local entanglement for a CDL tensor. 



square. In this paper, we do all our TEFR calculations 
using this scaling-SVD algorithm. 

It also clear that if the corner matrices in a CDL ten- 
sor have degenerate eigenvalues, the scaling-SVD cannot 
simplify such a CDL tensor. This is a desired result, since 
the degeneracy in the singular values may indecate the 
presence of topological order (see the discussion at the 
end of Appendix 4). 



2. Partition function and critical properties from 
the tensor network renormalization flow of the 
tensors 

In this section, we will give a more detailed discus- 
sion about the renormalization flow of the tensor. The 
more detailed discussion allows us to see the relation be- 
tween the fixed-point tensor and the energy spectrum for 
a quantum system. It also allows us to see the relation 
between the fixed-point tensor at the critical point and 
the scaling dimensions and central charge of the critical 
point. To be concrete, we will consider the renormaliza- 
tion transformation described by Fig. 2, which represents 
a single iteration step. Such an iteration step roughly 
correspond to two iteration steps described by Fig. 3. 

We start with a tensor network formed by No = N 
tensors = T. After i iterations, the tensor network 
is deformed into a tensor network of Ni T^-tensors. We 
formally denote the tensor network renormalization flow 
as 

T (i) _^ T (i+i)_ 

Those tensors satisfy the following defining property 

tTr^TW)^ =tTr(®T( i+1 )) JVi + 1 . 

The scales of the tensors T"W changes under each itera- 
tion. It is more convenient to describe the tensor network 
renormalization flow in terms of the normalized tensor. 
The normalized tensor is defined as 



(39) 
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FIG. 32: A square-lattice tensor network in a 2D space (or 
1+1D space time). 

that satisfies 

V r {i) - 1 

l,u 

In terms of the normalized tensor, after a tensor network 
rcnormalization transformation, a tensor network of Ni 
normalized T^-tensors is transformed into a tensor net- 
work of Ni + \ normalized T^ +1 ^-tensors: 

T« f l+1 T {l+1) ■ 

The N normalized T^-tensors and the N i+1 f i T { - i+r >- 
tensors have the same tensor trace 

tTr(®T®) Nt = /^\ +1 tTr(®T (j+1) ) Ar '+ 1 . (40) 

fi is a quantity that can be directly obtained from the 
numerical tensor network renormalization calculation. 

From (39) we see that the partition function at i th iter- 
ation can be expressed as a tensor trace of the normalized 
tensor 7"W up to an over all factor 

Zi = rf i tTr(®T«) jV *. (41) 

Note that by definition, the partition function is invariant 
under the tensor network renormalization transformation 
Zi = Zi + i. From (40), we find that Ti satisfies 

rf/^r^r (42) 

which allows us to calculate Ti from fi. We can keep 
performing the tensor network renormalization until Ni 
is of order 1. This way, we can calculate the partition 
function Z = Z t from (41) and (42). 

To calculate the density of the log partition function 
log(Z), we need to describe how tensor positions in the 
tensor network arc mapped into the physical coordinates 
r = [xi,Xi). Let X\ be the vector in physical space 
connecting the two neighboring tensors in the horizontal 
direction and X 2 the vector connecting the two neigh- 
boring tensors in the vertical direction (see Fig. 32). Let 
denote the coordinates of Xi as 

*«=(S£), a = 1,2. (43) 

The two by two matrix W describes the shape of the 
tensor network. 



We note that the tensor network is changed under each 
iteration. So W a p is also changed under each iteration. 
After the i th iteration, the shape of the resulting tensor 
network is described by wjfl, where W'f) = W a p is for 
the initial tensor network. Under the renormalization 
transformation described in Fig. 2, transforms in a 
simple way 

w (i+i) _ 2W { - €) . 

which implies Ai + \ = 4^4, . Under the renormalization 
transformation described in Fig. 4, transforms as 

W d+i) = W d) ^ ^ . 

The above recursion relations allow us to calculate WW 
for each step of iteration. Such information will be useful 
later. 

We see that the area occupied by each tensor TW after 
the i th iteration is 

A l = det(IU (i) ) 

where is the two by two matrix formed by W^p- It 
is also clear that 

AiN = A i+1 N l+1 = A tot 

is the total area of the system. 
Let us introduce 

z i = V\ /Ai . (44) 

From (42), we sec that Zj satisfies the following recursion 
relation 

z*i = ZifH? + \ (45) 

or 

log(z l ) = log(zo) + E^^. (46) 

3=1 3 

which will allows us to calculate the density of the log 
partition function 

A£ log(^) = log( Zi ) + A£ log (tTr(®T«)^) . 

Note that, after a number of iterations, we can reduce Ni 
to be of order 1. In this case log(zj) itself is the density 
of the log partition function in the large A to t limit. 

In the large iteration limit, and fi become inde- 
pendent of i. This allows us to introduce scale invariant 
tensor 

T (i) = 

inv / 

that satisfy 

rp(i) _^ rp(l+l) 
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in large i limit. We note that the invariant tensor has 
the following defining property 

From (40), we see that 7 must satisfy — j Ni + 1 / 4 JVl - +1 



/<+i 



or 



where we have used 



Ni=W i+1 . 



(47) 



(48) 



Thus the invariant tensor is simply given by 

= /r V3 ^ (i) (49) 

In the large i limit, /j approaches a constant and the 
recursion relation (45) implies that 

1 1 1 1 

\og(z i+1 f- + \ i+1 ) - logoff**) 

1 



3A, 



(Iog(/ i+ i) - log(/0) = oil/Ai). 



This result allows us to calculate the finite-size correc- 
tion to the partition function. For a tensor network 
of N tensors T, from the relation Z = tTr(<g)T) w = 
rf'tTr^TW)"*, wc find that the partition function is 
given by 

Z = zf«*tTr(®T®) Nt = {zif} *)^tTr(®7:^)^ 

When A tot and Ai are of the same order (ie when Ni is 
of order 1), we find that 

log(Z) = A tot log(z co ) + log (tTr(®T«)^) + o(l). 



time, the velocity is equal to 1 (or the critical system is 
isotropic in the 2D space). 

One way to obtain shape-dependent partition function 
(ie the A® ot part of the partition function) is to introduce 
four matrices associated with the invariant fixed-point 
tensor T- 



(0 



(M u %, d = £(T«) r „ rd , (M lr \ r = £(T«) r „,„, 

r u 
(M ldrU )ld,ru = (T^l)ruld, (M lurd )iu 7 rd = (T^l) ru ld- 

(50) 

The eigenvalues of the those matrices encode the infor- 
mation about that shape-dependent partition function. 

To understand the physical meaning of those eigenval- 
ues, let us consider a T^-tensor network that has a rect- 
angular shape Ni = L x M. The shape dependent part 
of the partition function is given by tTr((g)T i ^) 7Vi . To 

describe the shape quantitatively, we note that the T^- 
tensor network corresponds to a parallelogram spanned 

(i) 

by two vectors I = l a x a = LW\a x a and m = m a x a — 
MW^Xa, where x a , a = 1, 2 form an orthonormal basis 
of the space. The shape, by definition, is described by a 
complex number 

(j) _ mi + im 2 _ Wffi + iW 2 ff M 
h + ih ~ W$ + iW$ L ' 

Many critical properties of the critical point can be cal- 
culated from the r dependence of tTr(®T^) Ni . In par- 
ticular (according to conformal field theory) 71 



(tTr(®T«)^ 



f[('i"+'>'-H)Iii"- ( ' ) + i(hf-li')RcT (i) ] 



(51) 



n=0 



where z Q 



lim,- 



z ifi 



Also note that the invariant 



tensor T^l is independent of i in the large i limit (be- 
cause T-fnl approaches to a fixed-point tensor). We see 
that log(Z) contains a term A to tlog(z oc ) that is propor- 
tional to the system area A to t- log(Z) also contains a 

Atot independent term log ^tTr^T;^)^*^. The remain- 
ing o(l) term approaches to zero as A to t — * 00. 

It is interesting to note that although the 

log ^tTr((8)Ti nv ) Wi ^ term in the partition function 

does not depend on the total area A to t of the tensor 
network, it can depend on the shape of the network 
if the system is at the critical point. We can use this 
shape-dependent partition function to obtain the critical 
properties of the critical point. In the following, we will 
study the shape dependence of the partition function at 
a critical point. We will assume that in the 1+1D space 



where c is the central charge and {h^, h^} is the spec- 
trum of right and left scaling dimensions of the critical 
point. 

Since Tj^ is already a fixed-point tensor, we may eval- 
uate (51) by choosing L = 1: 

tTr(®T ; ( ^) M = Tr(M" d ) M 

_ ^2 e _2TM [( ft «+' l ^ra) Imr o i) + i ('i"-'^)Roi-o* ) ] 



n=0 



where 



(0 _ w£ + iWff 



(52) 



Using the above result, and some other results obtained 
with similar methods, we find that 
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>) + 


i(^ 


-^)Re(rf) 


(K +K- 


T2 )Im( 7o 


) + 






+ ^ - 












T2 )Im( ^ 


) + 




- ^)Re(^l) 



2tt 



2tt 
log(A^ d ) 

7T 

Iog(A«™) 



+ mod i, 
+ mod i, 



mod 2i, 



+ mod 2i, 



(53) 



where n = 0, 1, 2, ... and \% d , \%, X l ^ rd , and \ ldru are the 
eigenvalues of M ud , M lr , M lurd , and M w ™, respectively. 
Also 



Since /lo = 0, we see that the central charge c and the 
scaling dimensions /if , /if , ... can be calculated from the 

invariant fixed-point tensor T^. 

We like to stress that the invariant fixed-point ten- 
sor Ti nv can be calculated from (49) even away from the 
critical points. This allows us to calculate the central 
charge c and the scaling dimensions /if, /if , ... from Tj nv 
away from the critical point. For phases with short-range 
correlations, the resulting central charge c will be zero. 
However, the central charge will be non-zero when cor- 
relation length diverges (or energy gap vanishes). Thus 
plotting central charge as a function of temperature and 
other parameters is a good way to discover second or- 
der phase transitions. Also for phases with short-range 
correlations, a finite number of h^, are zero (which 
represent the degenerate ground states) and rest of them 
are infinite. 

We can also use the fixed-point tensors, to calculate the 
low energy spectrum of a quantum system, such as our 
spin-1 system. (29) From the tensor-network representa- 
tion of ID quantum Hamiltonian (see Fig. 20a), we see 
that the shape of the original tensor network is described 

by a VF-matrix W = . After i steps of itera- 

tion, the W-matrix becomes = 2 l W. From W^ l \ 

we find that is purely imaginary. When r}*^ oc i , the 
eigenvalues of the matrix — \og(M lurd ) will be propor- 
tional to the energy eigenvalues of our spin-1 system, as 
one can see from the expression of the partition function 
(51). From this result and using the fixed-point tensor, 
we can calculate the energy spectrum of quantum system. 



3. Phases and phase transitions from the 
fixed-point tensors 

In this section, we like to discuss how to define phases 
and phase transition more carefully and in a very general 
setting. For a system that depends on some parameters, 
for example (3 and V for our generalized loop-gas model, 
we can calculate the partition function Z((3, V) as a func- 
tion of those parameters. The phase and the phase transi- 
tions are determined from the partition function Z((3, V) 
in the following way: The singularity in the partition 
function Z(/3, V) marks the position of a phase transition. 
(This defines the phase transition.) The regions sepa- 
rated by the phase transition represents different phases. 
In other words, if a system V\) and another system 
(P2, V2) can be deformed into each other without phase 
transition, then the two systems are in the same phase. 
If there is no way to deform system (/?i,Vi) to system 
{02, V2) without encountering at least one phase transi- 
tion, then the two systems are in different phases. We 
like to point out that to avoid the complication of taking 
thermal dynamical limit, here we only consider systems 
with some translation symmetry. We also like to point 
out that the above definition is very general. It includes 
symmetry breaking phases, topological phases, as well as 
the phases with quantum order. 29 

However, the above definition is still incomplete. To 
define phases and phase transitions, it is also very im- 
portant to specify the conditions that the Hamiltonian 
or Lagrangian must satisfy. It is well known that with- 
out specify the symmetry of Hamiltonian, the symmetry 
breaking phases and symmetry breaking phase transi- 
tions cannot be defined. After we specify the conditions 
on the Hamiltonian, the deformations used in defining 
phases and phase transitions must satisfy those condi- 
tions. 

Since the conditions on the Hamiltonian or Lagrangian 
are part of the definition of phases and phase transitions, 
it is not surprising that the meaning of phases and phase 
transitions varies as we vary the conditions on the sys- 
tems. Let us illustrate this point using the generalized 
loop-gas model along the (3 — line. In fact, weather the 
fixed-point tensors T LL for V < V c and T Z2 for V > V v 
in the generalized loop-gas model describe two different 
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phases or not depend on the conditions on the system. 

Since we are using tensor network to describe our sys- 
tems, the conditions on the systems (or the Hamilto- 
nian/Lagrangian) will be conditions on the tensors. If 
the only condition on the tensors is that the tensors T ru di 
must be real, then T LL and T 7,2 describe the same phase. 
To show this, let us consider a family of tensors labeled 

T ruld = RrM9) R uu>(9) T r> 2 u>l>d> R l'l(Q) R d>d(8) 

where R(0) = " „ . We see that the tensors 

v ' \smv — cos 9 J 

T e connect T Z2 = T e=0 and T LL = T 9=7r / 4 as we change 
9 from to 7r/4. Since R{9) is a field redefinition trans- 
formation, the partition function for T e does not depends 
on 9 and is hence a smooth function of 9. Therefore, T LL 
and T 7,2 describe the same phase. 

On the other hand, if we specify that the tensors must 
satisfy the closed loop condition (23) (or (27)) and the 
Z 2 symmetry condition (26), then T LL and T 7,2 describe 
two different phases. The continuous deformation T e can 
no longer be used. This is because T does not satisfy 
(23) if < 9 < tt/4. 

Before ending this section, we like to point out that due 
to field redefinition ambiguity in the TEFR calculation, 
some time, a single phase may correspond to a group 
fixed-point tensors that are related by field redefinitions. 
To understand this point more concretely, let us consider 
in more detail the V > V c phase described by the fixed- 
point tensor T Z2 . In the decomposition of rank-4 tensor 
T into two rank-3 tensors S 2 and S4: 

Truld = y^XSz)ldk(S4)ruk 

a 

we may have a ambiguity that T can also be composed 
into two different rank-3 tensors S' 2 and S' 4 

Truld = y ^X ^2)ldk{S^)ruk 

a 

where S 2 and S 2 , S'4 and S' 4 are related by an invertible 
transformations 

(S 2 )ldk = Akk' (S' 2 )ldk' {Si)ldk = {S' A )ldk'{A 1 )k'k- 

Here A represents the field redefinition ambiguity in the 
TEFR renormalization approach. 

If the tensors T and S must satisfy certain conditions, 
such as the closed loop condition (23), those conditions 
will limit the ambiguity A. For the V > V c phase, the 
fixed-point tensor is T Z2 . The closed loop condition (23) 
limits the ambiguity of the decomposition of T 7,2 to a 
diagonal A matrix. The field redefinition (19) does not 
change T 7,2 when A and B in (19) are diagonal matrices. 
This way, we show that the fixed-point tensor in the V > 
V c phase has no field redefinition ambiguity. Similarly, 
for the V < V c phase, the fixed-point tensor T LL also 
has no field redefinition ambiguity, if we implement the 
closed loop condition (23). Those results are confirmed 
by our numerical calculation. 



4. A detailed discussion of the renormalization flow 
of tensors 

Let us describe the TEFR flow in our spin-1 model in 
more detail which allows us to see how an initial time- 
reversal and parity symmetric tensor T flows to the var- 
ious fixed-point tensors. Let us first consider only the 
SVDTRG transformation described in Fig. 4. 

It is most convenient to describe the tensor in term of 
its singular values X s under the SVD decomposition (see 
(3)). At first a few steps of iterations, the singular values 
A s are in general not degenerate. If the initial tensor T 
is in the Haldane phase, then after several iterations, the 
four largest singular values become more and more de- 
generate. If the initial tensor T is in the TRI phase, then 
after several iterations, the largest singular values remain 
non degenerate. We see that the Haldane phase and the 
TRI phase can be distinguished from their different ten- 
sor flowing behaviors under SVDTRG transformation. 

Now let us consider the tensor flow under the TEFR 
transformation in Fig. 7. Again at first a few steps of 
iterations, the singular values A s are in general not de- 
generate. Also the tensors in first a few steps are quite 
different from CDL tensor. As a result, the scaling-SVD 
operation in Fig. 7(b)— >(c) is a null operation which does 
not simplify the tensors. After several iterations, the ten- 
sors become closer and closer to CDL tensors. So after 
a certain step of iterations, the scaling-SVD operation 
becomes effective which further simplify the tensor. 

If the initial tensor T is in the TRI phase, even after the 
tensor T becomes CDL tensor, the singular values X s are 
still not degenerate. This implies that the corner matri- 
ces Mi in the CDL tensors (see (6)) have non-degenerate 
eigenvalues. The entanglement filtering operation in Fig. 
7(b)— >(c) actually simplify the tensors by removing the 
smaller eigenvalues in the corner matrices Mi. So the 
scaling-SVD operation reduces the corner matrices Mi 
to a trivial one (Mj),j/ = Si^Si^. This way the scaling- 
SVD operation reduces the initial tensor to the trivial 
dimension-one tensor T™. 

If the initial tensor T is in the Haldane phase, after the 
tensor T becomes CDL tensor after a few steps of TEFR 
iteration, the largest four singular values X s also become 
degenerate. This implies that the largest two eigenval- 
ues of the corner matrices Mj in the CDL tensors are 
also degenerate. In this case, the scaling-SVD opera- 
tion cannot reduce such a degenerate corner matrices to 
a trivial one. But it can still simplify the corner matri- 
ces to a certain degree by removing other smaller eigen- 
values. This way, the combinations of the SVD coarse 
graining in 7(d)— >(e) and the entanglement filtering in 
Fig. 7(b)— >(c) reduce the initial tensor T to a simple 
fixed-point tensor T Haldanc . Therefore the Haldane phase 
and the TRI phase can be distinguished by their differ- 
ent fixed-point tensors under the TEFR renormalization 
transformation. 

We like to point out the emergence of the 4- fold degen- 
erate singular values requires the initial tensors to satisfy 
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the time-reversal and parity symmetries. If we start with 
an initial tensor T that has time-reversal symmetry but 
not parity symmetry, our numerical results show that the 
largest singular values do not become degenerate after 
many iterations. In this case the entanglement filtering 
operation through the scaling-SVD always reduces the 
initial tensor T to T TRI . This implies that, without the 
parity symmetry, the Haldane phase is the same as the 
trivial phase TRI. 

We see that the emergence of 4-fold degenerate sin- 
gular values is closely related to the Haldane phase - a 
topologically ordered phase. In Ref. 72, an entanglement 
spectrum is introduced to identify topological order. We 
believe that the emergence of 4-fold degenerate singular 
values is closely related to the degeneracy in the entan- 
glement spectrum. 

We also like to point out that, only away from criti- 
cal points, the tensors flow to CDL tensor. Near critical 
points, tensors are always quite different from CDL ten- 
sors. As a result, the entanglement filtering operation in 
Fig. 7(b)— >(c) is a null operation which does not simplify 



the tensors. So near a critical point, the SVDTRG in Fig. 
4 and the TEFR in Fig. 7 produce the same result. 



5. Projective realization of symmetries on 
fixed-point tensors 

We like to point out that in the original tensor net- 
work of T for the spin-1 system (30), the time reversal 
and parity symmetry are given in (32) and (33). How- 
ever, as discussed near the end of the Appendix 3, the 
TEFR transformation has a field-redefinition ambiguity. 
As a result, the time-reversal and parity transformations 
contain a field-redefinition ambiguity, when act on the 
coarse grained tensor network. Similarly, the fixed-point 
tensor also contain such field-redefinition ambiguity. If 
we choose a basis that the fixed-point tensor in the Hal- 
dane phase is given by T Haldano = T(a y , a y , a y , a y ), the 
time reversal and parity transformations that act on the 
fixed-point tensor are given by 



Parity: T -> T p , where {T p ) ruU = ^ W rrl W uu ,W w W ddl T u , r , d r l 



r'u'l'd' 



Time reversal: T —>T t , where (T t ) ruW = T dlv 



r 



(55) 
(56) 



Here W is a field-redefinition transformation 



W 



(\ 0\ 

10 

10 

\0 1/ 



(57) 



The time-reversal and parity transformations on the 
coarse grained tensor network and those on the ini- 
tial tensor network (see (32) and (33)) differ by a 
field-redefinition transformation. The fixed-point tensor 
^Haldane j g mvar j an t under the above projective symme- 
try group (PSG) transformations. 29 We also numerically 
confirm that the fixed-point tensor T Haldanc is stable if it 
is perturbed by a tensor ST that is invariant under the 
same PSG defined through (55) and (56). We note that 
the PSG is used to define quantum order (a generaliza- 
tion of topological order) in Ref. 29. The tensor network 
renormalization naturally allows the generalization from 
topological order to quantum order. 

We like to remark that the field-redefinition trans- 
formation should be derived by keeping track of the 
parity and time reversal symmetry at each TEFR flow 
step. As a result, the field-redefinition transforma- 
tion W should satisfy (55) for the fixed point tensor 



yHalda 



T(a v ,a y ,a y ,a y ). However, not all the W 



that satisfy (55) for the fixed point tensor T Halda 
T(a y 7 a y ,a y ,<j v ) correspond to the parity symmetry for 



the original tensor (33). For example: 



w 



/0 IN 
0-100 

0-10 

\1 0/ 



(58) 



also satisfies Eq. (55), but this is not the correct 
field-redefinition transformation for parity symmetry of 
T Haidane = T{a y , a y , a y , a y ). In particular, an perturba- 
tions with the wrong parity symmetry: 

ST ruld = J2 W' rr ,W' uu ,Wl v W dd ,5T ul 

r'u'l'd' 



r'd'V 



(59) 



and time reversal symmetry (56) can destabilize the 
fixed-point tensor T rialdano = T(a y , cr y , cr y , a y ). 

For the fixed point tensor of Haldane phase, the field- 
redefinition transformations that correspond to parity 
and time reversal symmetries are also base dependent. 
If way choose T Haldano = T(a°, cr°, a , <7°), then the 
correct field-redefinition transformation corresponding to 
the parity symmetry should be W not W. We can also 
choose other bases that T Haldane is invariant in a trivial 
way under parity transformation but invariant up to a 
field-redefinition under time reversal transformation. 
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